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Abstract. - We numerically analyze Non-Oberbeck-Boussinesq (NOB) effects in two-dimensional 
Rayleigh-Benard flow in glycerol, which shows a dramatic change in the viscosity with temper- 
ature. The results are presented both as functions of the Rayleigh number Ra up to 10 s (for 
fixed temperature difference A between the top and bottom plates) and as functions of A ( "non- 
Oberbeck-Boussinesqness" or "NOBness") up to 50K (for fixed Ra). For this large NOBness the 
center temperature T c is more than 5K larger than the arithmetic mean temperature T m between 
top and bottom plate and only weakly depends on Ra. To physically account for the NOB devia- 
tions of the Nusselt numbers from its Oberbeck-Boussinesq values, we apply the decomposition of 
Nunob/Nuob into the product of two effects, namely first the change in the sum of the top and 
bottom thermal BL thicknesses, and second the shift of the center temperature T c as compared 
to T m . While for water the origin of the TVm deviation is totally dominated by the second effect 
(cf. Ahlers et al. J. Fluid Mech. 569, 409 (2006)) for glycerol the first effect is dominating, in 
spite of the large increase of T c as compared to T m . 
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Introduction. — In most theoretical and numer- 
ical studies on Rayleigh-Benard (RB) convection, the 
Oberbeck-Boussinesq (OB) approximation [1, 2] is em- 
ployed, i.e., the fluid material properties are assumed to 
be independent of temperature T except for the density in 
the buoyancy term which is taken to be linear in T. The 
problem has two control parameters, namely the Rayleigh 
number Ra = f3gL 3 A/(kv) (here (3 is the thermal ex- 
pansion coefficient, g the gravitational acceleration, L the 
height, A the temperature difference between bottom and 
top plates, k the thermal diffusivity, and v the kinematic 
viscosity), and the Prandtl number Pr = v/k. For the OB 
case the mean temperature profile shows top-bottom sym- 
metry. However, in real fluids, if A is large, this symmetry 
no longer holds due to the temperature dependences of 
the material properties. Thus, for given fluid, A appears 
as an additional control parameter, which characterizes 
the deviations from OB conditions, leading to so called 
Non-Oberbeck-Boussinesq (NOB) effects. The NOB sig- 
natures can be quantified by (i) a shift T c — T m of the 
bulk (or center) temperature T c from the arithmetic mean 



temperature T m between the bottom and top plates) and 
(ii) by the ratio of the Nusselt numbers Nu^o b /Nuo b in 
the NOB and OB cases, which deviates from one. Both 
quantities have been measured in the large Ra regime for 
helium [3], glycerol [4], ethane [5], and water [6] as func- 
tions of the NOB-ness A. 

As shown in Ahlers et al. [6] the Nusselt number ratio 
Nunob/Nuob can be connected to T c by the identity 
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Here the labels on material properties indicate the tem- 
perature at which they arc taken, e.g. n t = n(T t ) etc. 
At = T c — T t and Af, = TJ, — T c denote the temperature 
drops over the top and bottom thermal boundary layers, 
and Xf and Xf indicate their thicknesses, based on the 
temperature slopes at the top and bottom plates, respec- 
tively. Xq B is the thermal BL thickness in the OB case, 
both at top and at bottom. The factor Fa can be calcu- 
lated from the temperature dependences of the material 
properties immediately, once T c is known. Remarkably, 
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Ahlers et al. [6] experimentally found that for water 

F X « 1. (2) 

This has been confirmed by numerical simulations of 2D 
NOB Rayleigh-Benard convection in ref. [7] . If the relation 
eq.(2) holds, the Nusselt number ratio N nu°b a l reaa -y 
lows from the center temperature T c , which for water can 
be calculated within a generalized boundary layer theory 
introduced in ref. [6]. 

The objective of this letter is to answer the appar- 
ently important question whether the relation Fj w 1, 
meaning that the sum of the boundary layer thicknesses 
stays the same as in the OB case also under NOB con- 
ditions, Xf + Xf w 2Aq B , is more generally valid, i.e., 
if it holds for other liquids too. We therefore have per- 
formed (two-dimensional) NOB simulations with glycerol 
as the working fluid. For glycerol the kinematic viscosity 
dramatically depends on temperature, i.e., one should ex- 
pect large changes of the boundary layer thicknesses at top 
and bottom. For instance, v decreases from 1759 mm 2 /s 
to 52.5 mm 2 /s if the temperature increases from 15°C to 
65°C. Another advantage of considering glycerol is the ex- 
istence of experimental data for the center temperature 
(see ref. [4] ) for comparison (but not for the Nusselt num- 
ber modification) . Our main result will be that the sum of 
the boundary layer widths is indeed changed under NOB 
conditions, i.e., relation (2) does not hold for glycerol. Its 
validity for water thus turns out to be coincidental, due to 
the specific temperature dependences of its material pa- 
rameters for the chosen temperatures in the experiments 
of ref. [6] . 

Note that the fluid flow in glycerol is very different from 
that in water at the same Ra. Due to glycerol's huge 
Prandtl number of about Pr m 2500, the transition range 
between the onset of convection at Ra c and the loss of 
spatial coherence in the flow is much more extended than 
for water or air, whose Prandtl numbers are of order one. 
While in air and in water this transitions range extends 
to about Ra w 5 ■ 10 7 to 10 8 , only beyond which there 
is turbulent convection in the bulk, this range extends to 
much larger Ra in glycerol, namely to Raylcigh numbers 
of order Ra m 10 12 . Since the numerical calculations cover 
the range up to Ra « 10 s only, all results refer to a fluid 
flow still having coherent structures. 

To quantify these statements we use an averaged Kol- 
mogorov length tjk as a measure for the scale of coherent 
structures in the flow, more precisely 

tcoh = 10 m = 10 {^/e u ) 1/A . (3) 

Here e u is the volume average of the energy dissipation 
rate of the flow for which the well-known exact relation 
e u = v z L^ i Pr~ 2 Ra{Nu — 1) holds. With this we obtain 

l coh /L = 10 Pr 1 ' 2 (Ra(Nu - 1))" 1/4 (4) 

as an estimate for a volume-averaged relative coherence 
length. Taking Nu(Ra, Pr) from the unified theory of 
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Fig. 1: (color online) The coherence length £ co h in multiples 
of the cell size L versus Ra number for three fluids. For air 
and water it is of order 0.1 near Ra « 10 7 to Ra « 10 s . For 
glycerol this is reached much later; turbulent heat convection is 
expected beyond Ra ~ 10 12 only. Gray shaded region indicates 
the developed turbulent regime. Lines are derived from the 
unified theory of ref. [8], symbols correspond to the Rayleigh 
numbers of the present OB numerical simulations. 

refs. [8], one thus obtains an estimate of the coherence 
length as a function of Ra and Pr from eq. (4), see figure 
1. The main features of the coherence length are (i) its 
pronounced explicit dependence on Pr (the implicit de- 
pendence via Nu is only weak). It is by about a factor 
\/2500 = 50 larger for glycerol than for gases or water, 
(ii) Its i?a-depcndcnce is approximately £ co h cx Ra~ 3 . 

General description of numerical simulation. 

We numerically solve the incompressible (diUi = 0) 
Navicr-Stokcs equations 

p m (d t Ui + UjdjUi) =-dip + djir/idjUi + diUj)) 

+ g{pm - p)S i3 , 

and the heat-transfer equation 

(d t T + u 3 d 3 T) = dj{hBjT). (6) 

The temperature dependence of the dynamic viscosity 
T)(T), the heat conductivity A(T), and the density p are 
experimentally known for glycerol. They are given in the 
appendix of reference [6]. As justified in that reference, 
we can assume the isobaric specific heat capacity c p and 
the density p in the time derivatives of the material pa- 
rameters to be constant at their values p m and c Pj7Ti at the 
arithmetic mean temperature T m . We vary the Raylcigh 
number Ra up to 10 8 and the level of the NOBness A up 
to 50K. 

The container is two-dimensional (2D, no y- 
dependence), has height L, and aspect ratio 1. The 
flow is wall-bounded, i.e., we use no-slip boundary condi- 
tions at all solid boundaries: Ui = at the top (z = L) 
and bottom (z = 0) plates as well as on the side walls 
x = and x — L. For the temperature at the side walls 
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heat-insulating conditions are employed and Tf, — T t = A 
is the temperature drop across the whole cell. The 
Raylcigh number is defined with the material parameters 
taken at the mean temperature T m , i.e., Ra = /9 " g ^ A . 
The arithmetic mean temperature is fixed at T m = 40°C. 
We vary the Rayleigh number by varying the height L of 
the box, while the NOBness is changed by varying the 
temperature drop A. Note that in the buoyancy term 
in eq.(5) the full temperature dependence of the density 
is taken into account, rather than employing the linear 
approximation p(T) — p m = p m (3{T — T m ) only. (Never- 
theless, the Raylcigh number is defined as usual with the 
linear expansion coefficient of the density with respect 
to temperature, taken at T m , namely /3 m = — ~5rlT m -) 
The Prandtl number is defined as Pr = v m jK m \ for 
glycerol at the chosen temperature T m its value is 
Pr = 2495. The basic equations are directly solved on 
the two-dimensional domain by means of the fourth-order 
finite difference method. For a detailed description of the 
simulation method as well as its validations, see ref. [7]. 

One may worry if two-dimensional simulations are suf- 
ficient to reflect the dynamics of the three-dimensional 
RB convection. For convection under OB conditions this 
point has been analyzed in detail in ref. [9] and earlier 
in refs. [10-13]. The conclusion is that for Pr > 1 var- 
ious properties observed in numerical 3D convection and 
in experiment are well reflected in 2D simulations. This 
in particular holds for the BL profiles and for the Nussclt 
number. Since the focus of this paper is on the difference 
between OB and NOB convection, the restriction to 2D 
simulations seems to be even less severe, as NOB devia- 
tions are expected to be similar in both 2D and 3D simula- 
tions and remaining differences to cancel out in quantities 
such as T c — T m or Nunob/Nuob- 

Results and discussions. 

Large scale flow dynamics and temperature snapshots. 
In the steady flow regime (Ra < 1.5 • 10 5 ) a single large- 
scale circulation role develops, which however disappears 
in the unsteady flow regime (Ra > 1.5 • 10 5 ) and does not 
reappear up to the largest accessible value Ra = 10 8 of 
the present study. Even if we start the simulation with an 
artificial single roll, the large-scale circulation disappears 
in the course of time and then isolated plumes (as shown 
in figure 2) dominate the flow. This feature holds for both 
cases, OB and NOB, and is qualitatively different from the 
observations in 2D (OB and NOB) simulations in water 
(see ref. [7]). We attribute this to the much larger spatial 
correlations in glycerol as addressed above. Note that in 
experiment (ref. [4]) for larger Ra = 2.3 ■ 10 s a large- 
scale 3D circulation role has been observed for glycerol. 
The different behavior between the present DNS and the 
experiment could either be due to the smaller Ra or to the 
two-dimensionality in the simulation. 

Typical temperature snapshots are shown in Figure 2. 
As observed in experiments, refs. [4,6], the NOB convec- 
tion is characterized by an enhancement of the bulk tem- 



1.0 



0.8 



0.0 



'wmh^ j 



mil' 1 ''' 
o.6 y\\f"~--''iiii\\\w 

0.4 p ~"«wfc\\\\v-- 
'iin 

0.2 




0.0 



0.2 



0.4 0.6 
x/L 



0.8 



1.0 





20 



40 



60 o 



c 



Fig. 2: (color online) Snapshots of the velocity and temper- 
ature fields for Ra = 10 s at T m = 40° C. The upper panel 
corresponds to the OB case (T-independent material param- 
eters), the lower one to the NOB case, both with A = 40K. 
The temperature color scheme is the same in both cases. In 
the NOB case a strong temperature enhancement of the center 
is clearly visible. 



perature T c , and a top-bottom asymmetry of the thermal 
BL thicknesses. Due to the large variation of the glycerol 
viscosity (the viscosity ratio reaches as much as vt/vb ~ 16 
at A = 40K), the more viscous cold plumes from the top 
BL arc much less mobile than the warmer plumes from 
the bottom BL. This results in a significant increase of T c 
as compared to the water case. 

Mean temperature profiles and center temperature. To 
quantify the enhancement of the bulk temperature T c , the 
temperature profiles for Ra = 10 s are shown in Fig. 3. 
Again, a strong asymmetry between top and bottom is ob- 
served: Due to the more mobile bottom plumes the center 
temperature T c is significantly larger than T m . 
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Fig. 3: (color online) Mean temperature profiles for glycerol at 
Ra = 10 s in the OB case (dashed) and in the NOB case with 
A = 40A' (solid). (In both cases T m = 40°C, same Ra, same 
A, but T-independent (OB) or T-dependent (NOB) material 
parameters, respectively.) In the NOB case the strong temper- 
ature enhancement of the center temperature T c by about 5 K 
becomes visible (relative shift ss 0.12). 



To demonstrate this the center temperature shift T c — 
T m (normalized by A) as function of the Raylcigh number 
Ra and of the NOBness A is shown in Figure 4. Except 
for small Raylcigh numbers just above onset of convection 
and in a region around Ra w 2 • 10 5 just above the onset of 
unsteady motion, the bulk temperature shift (T c — T m )/A 
is rather independent of Ra. The tiny increase between 
Ra = 10 7 and 10 s however is beyond the statistical error- 
bars. For comparison, the prediction of the NOB BL the- 
ory given in ref. [6] and the shift for the non-convective 
state (i.e., purely conductive heat transport, driven by the 
temperature gradient only) are shown. Though the NOB 
BL theory from ref. [6] is not applicable here due to the 
lack of a large scale wind, it gives the correct qualitative 
trend for the shift (T c — T m )/A. We also included ex- 
perimental data measured at Ra = 2.3 • 10 s (taken from 
ref. [4]) in an aspect ratio 1 cylindrical container. Though 
for that case a large scale convection role has been ob- 
served, the agreement with the 2D numerical simulations 
is reasonable. 

Nusselt number. The key question on NOB effects is: 
How do they affect the heat flux, i.e., the Nusselt num- 
ber? For water we could address this question within an 
extended BL theory, cf . ref. [6] , but only thanks to the ex- 
act relation eq.(l) and the experimental input F\ ps 1, see 
relation (2), because then only Fa is needed to calculate 
the NOB deviations in the Nusselt number ratio, and Fa 
is accessible within the extended BL theory, since it fol- 
lows directly from T c . But here, with glycerol as working- 
fluid, wc find that F\ fa 1 does not hold, as demonstrated 
in Fig. 5. In contrast to water, for glycerol the main A- 
dcpcndcncc of Nunob/Nuob — F\ ■ Fa is due to the 
A-depcndence of Fx while the factor Fa is basically 1 for 
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Fig. 4: (color online) Relative deviation (T c — T m )/A of the 
center temperature T c from the arithmetic mean temperature 
Tm for glycerol versus Ra at fixed A = 40K (upper) and versus 
the NOBness A at fixed T m = 40°C and various values of Ra 
(lower). The experimental data points (denoted by x) are 
measured at Ra — 2.3 • 10 8 and were taken from ref. [4]. For 
comparison the T c shift obtained from BL theory (upper solid 
lines) and for the case of no convection (lower solid lines) are 
also plotted. 



all A. This qualitative difference between glycerol and wa- 
ter in the origin of the Nusselt number modification also 
means that the experimental finding F\ ~ 1 for water at 
T m = 40F and Ra in the range of 10 8 - 10 10 , see ref. [6], 
is merely accidental and not a general feature of the RB 
flow under NOB conditions. 

Both NOB responses, the shift of the center temper- 
ature T c and thus Aj, ^ A t as well as the shift of the 
BL thicknesses Ajj t , are determined by the full nonlinear 
dynamics, in glycerol as well as in water. The T c -shift in 
glycerol is even larger 6.5K) than in water IK). The 
same is expected for the A^ t -shifts. But the differences in 
the temperature drops A;,^ enter via Fa; here they are 
weighed with the explicit temperature dependence of the 
material parameter n(T). Since the thermal diffusivity 
changes only minutely in glycerol, K(,.t/K m — 1 w ±0.01, 
the factor Fa stays near Fa ~ 1 despite the large T c re- 
sponse, cf. Figs. 5,6. This does not happen in Fx; here 
the full changes of A^ t enter. Because of the very strong 
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Fig. 5: (color online) Nusselt number ratio Nunob /Nuob — 
F\ ■ Fa together with its contributing factors Fx and Fa versus 
A for fixed Rayleigh numbers, (a) Ra — 10 4 , (b) Ra = 10 7 , 
and (c) Ra — 10 8 . As always, the working liquid is glycerol at 
Tm = 40° C. The dashed lines correspond to Fa resulting from 
the NOB BL theory of ref. [6]. 



and in particular nonlinear temperature dependence of v 
the thicknesses of the BLs change significantly and also 
quite differently in magnitude at the bottom and the top 
BLs, because \J v t ~jvh ~ 2.4 due to the strong nonlinear 
T-dependcncc of v{T). Therefore the sum Xf + Xf no 
longer is equal to 2Aq B . For water, instead, the domi- 



nantly linear A^-NOB modifications are opposite in sign 
and nearly cancel in the sum of the NOB thicknesses, giv- 



ing Xr + Af 



2Aq B or F\ 



1. Thus in glycerol we 



have Fa ~ 1 and the Nu changes are dominated by Fa, 
while in water it is Fa ~ 1 and the NOB effects in Nu 
are determined dominantly by Fa (which is given by the 
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Fig. 6: (color online) The Nusselt number ratio 
Nunob /Nuob — Fx ■ Fa (upper) and the constituting 
factors Fa and Fa individually (lower) versus Ra for fixed 
NOBness A = 40A" (glycerol at T m = 40°C). Note the 
dramatic NOB effect at Ra w 2 • 10 5 ; this happens still in 
the pattern forming range, far below the turbulent high Ra 
region. We are not aware of its experimental verification. 



temperature shift alone). 

Figure 6 shows that the dependences of Nunob/Nuob 
and F\ on the Rayleigh number Ra are non-monotonous. 
We consider this as due to the nontrivial evolution of var- 
ious coherent flow patterns with increasing Ra. In par- 
ticular, as shown in figure 5, for Ra = 10 the function 
Fa (A) shows a qualitatively opposite behavior to that for 
water, namely Fa increases with increasing A and reaches 
as large a value as 1.017 at A = 50K. A consequence of 
our finding is that in general Nunob/Nuob = Fx ■ Fa 
cannot be calculated within the extended BL theory in- 
troduced in ref. [6] , even if a large-scale wind has formed: 
Within BL theory only the factor Fa can be calculated 
but not the factor Fa, for which in general one cannot 
assume Fx ~ 1. 

We finally present our results for the Nu number itself 
as a function of Ra, see Fig. 7, both for the OB and the 
NOB case. The inset shows the local scaling exponents. 
When applying the unifying theory of refs. [8], it is 0.306 
at Ra — 10 s and Pr = 2500, consistent with our numerical 
findings. This local slope practically does not change in 
the NOB case. 
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Fig. 7: (color online) The Nusselt number Nu for glycerol ver- 
sus Ra under OB (dashed line) and NOB (solid line) condi- 
tions. In both cases T m = 40°C and A = 40K. OB is pro- 
vided by keeping the material parameters artificially constant 
with T. We have also included the available data from ref. [4]. 
Logarithmic slope dlog(iVii)/d log(i?a) is plotted in the inset 
and the line corresponding to the exponent 0.297 measured in 
ref. [4] is also shown. 



Summary and conclusions. — In summary, for glyc- 
erol both the center temperature T c and the Nusselt num- 
ber Nu of the 2D numerical simulations are in good agree- 
ment with the available experimental data of ref. [4] . The 
experimental finding by Ahlers et al. [6] of a "thermal-BL- 



thickness sum rule" for water, F\ 



1 or X s b l 



-a; 



is shown to be incidental and seems due to the specific tem- 
perature dependence of the material parameters of water 
at 40°C. Apparently this cannot be generalized to other 
fluids (or other mean temperatures), as our analysis of 
RB convection in glycerol has shown. While for water the 
Nusselt number modification Nu nob /Nuo b is due to the 
modified temperature drops over the BLs, represented by 
_Fa, as shown in refs. [6,7], for glycerol it is governed by 
the variation of the BL thicknesses, namely by F\. This 
can be attributed to the strong and nonlinear temperature 
dependence of v{T). 
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